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The magnetic-field dependence of the energy spectrum of GaAs doped with nitrogen impurities 
is investigated. Our theoretical model is based on the phenomenological Band Anticrossing Model 
(BAC) which we extended in order to include magnetic field and electron - phonon interaction. Due 
to the highly localized nature of the nitrogen state, we find that the energy levels are very different 
from those of pure GaAs. The polaron correction results in a lower cyclotron resonance energy as 
compared to pure GaAs. The magneto-absorption spectrum exhibits series of asymmetric peaks 
close to the cyclotron energy luo c . 
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I. INTRODUCTION 

Contemporary epitaxial growth techniques provide the 
possibility for synthesis of high quality semiconductor al- 
loys and/or elemental materials. Recently, there has been 
growing interest in N doped GaAs, due to its possible ap- 
plications for long wavelength optoelectronic devices [l|. 
In GaAsi-zNu substitutional nitrogen has the same va- 
lence state as As, thus forming an isoelectronic impurity. 
Experimental data have unambiguously shown that the 
effects of nitrogen incorporation are at least three-fold: a) 
reduction of the fundamental band- gap 0-(j| i b) change 
in the electron effective mass [a, Q; and c) decrease in the 
electron mobility Q. 

Substitutional nitrogen in GaAs forms a resonant level, 
above the conduction minimum[8] (El = 0.23eV), hav- 
ing Ai (spherical) symmetry. Since N(2s 2 2p 3 ) has the 
same valence state as As(4s 2 4p 3 ) they differ mainly in 
their local pseudo-potentials resulting in interaction that 
is predominantly short-range [9, 10j. It is found that the 
position of this resonant level, El, does not change for 
nitrogen concentration up to x — 3%, and the bowing 
of the conduction minimum follows a simple square root 
like law for small x. Therefore one may expect that a 
phenomenological approach using perturbation theory is 
sufficient to explain, for instance, the reduction of the 
band-gap. However, this is true only in part, as the sig- 
nificant difference in bond lengths between Ga-As and 
Ga-N affects next-nearest neighbours [Tl|, so that the 
impurity potential has a part which is of intermediate 
range. What is more, for higher concentration of nitro- 
gen (x > 1%), GaAs:N is classified as a semiconductor 
alloy, when the use of Virtual Crystal Approximation 



(VCA) is questionable, and partial colla pse of the Bril- 
louin zone is expected. Some authors [12, LL3| distinguish 
three ranges for the molar concentration: 1) ultradilute 
(x < 0.01%), 2) dilute (0.01% < x < 1%), 3) semicon- 
ductor alloy (x > 1%). In this paper we will confine our 
treatment of GaAs doped with moderately low concentra- 
tion of N (a; = 0.08%) which corresponds to case 2). We 
will investigate the effect of a high magnetic field and of 
electron-phonon interaction on the energy spectrum and 
on the magneto- absorption spectrum of GaAsi^N^ that 
has not been considered up to now. 

The paper is organized as follows. In Sect. [TT| we 
present our theoretical formalism which is based on the 
phenomenological Band Anticrossing Model (BAC). It is 
shown how to derive the energy spectrum in the pres- 
ence of magnetic field in bulk GaAs:N. In the following 
section, Sect.|TTT]we discuss the influence of the electron- 
phonon interaction on the fundamental transition energy, 
i.e. the difference between the first two Landau levels. 
In Sect. IIV1 theoretical estimates of the absorption spec- 
trum are given for two different values of temperature 
and magnetic field. In the last section, Sect. [V] we sum- 
marize the results and present our conclusions. 



II. THEORETICAL FORMALISM 

The Band Anticrossing Model (BAC) will be employed 
to determine the energy levels in GaAs:N under an ap- 
plied magnetic field B with inclusion of electron-phonon 
interaction. Within this model the interaction between 
the nitrogen state and the conduction states is charac- 
terized by a single value Cn, to be defined later. This 



is possible in dilute GaAs:N where the overlap of neigh- 
boring nitrogen wavefunctions is negligible. Then the to- 
tal wavefunction can be written as a linear combination 
of extended states ^ c of pure GaAs, and the impurity 



wavefunction ^ 
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It is convenient to represent the conduction band wave- 
function \£>q(t) in the basis of Wannier wavefunctions, 
ac (r — Rj ) , centered around the sites of the crystal Rj 



*W(r) = -^]To c (r-R,K k ^ 



(2) 



In the above expression M is the number of Ga-As 
pairs in arbitrary large volume of the crystal. The matrix 
element between the nitrogen state and conduction states 
of pure GaAs is 



n will refer to pure Landau levels, while n\ and n,2 will 
refer to lower and upper subband branches (n— and n+ 
in Eq. ([7])). In the absence of impurities, the conduc- 
tion (extended) wavefunctions have the form of a linear 
harmonic oscillator, and for the magnetic field orienta- 
tion along the z axis and choosing the Landau gauge 
A = (—By, 0, 0), their explicit form is as follows 
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Notation H n corresponds to the n th order Hermitian 
polynomial, and yo is the y— coordinate of the cen- 
ter of the orbit, while the normalization constant de- 
pends on the cyclotron orbit l c , and is given by N n — 

(V^ c 2"n!)- 1/2 . 

The final wavefunction for the n[ h lower state has the 
form: 



tt ni (r)=a ni tf^(r)+/3 ni tf L (r). 



(9) 
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where x is the concentration of nitrogen impurities in the 
semiconductor. Eqs. ([I])-© give rise to a two-level like 
secular equation [9J, 



E-E c (k) ^x~C N 
ViCW E - E L 



= 0, 



(4) 



where Ec{k) is the energy dispersion of the conduction 
band, and El is the N-impurity level. The solutions of 
Eq. (0} are 



E± = - { E c + E L d 



(E C - E L f 



4xC 2 N 
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The value of Cn is estimated to be around 2.7eV by 
fitting to experimental data 9]. The envelope wave func- 
tion of a nitrogen impurity can be approximated by a 
Gaussian-like function having only one parameter, the 
localization radius a, 
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In the presence of a magnetic field B, the conduction 
band splits in a series of Landau-like levels whose energy 
spectrum is modified by the presence of the short range 
impurity potential AV A (r), 
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FIG. 1: Energy levels in GaAs:N with nitrogen concentration 
of x — 0.08%, as a function of magnetic field B, neglecting the 
interaction with phonons. The levels are derived from Landau 
levels, whose indices are shown near the different curves. It 
is obvious that the spectrum consists of two branches, one of 
which is below the nitrogen level, El, while the another one 
is above. 

Fig. [1] shows the energy levels in GaAs:N with nitro- 
gen content of x = 0.08%, vs magnetic field strength. 
The interaction with optical phonons is neglected for the 
moment to emphasize the influence of the magnetic field 
alone. The spectrum splits into two parts with respect to 
the isolated, localized level El- The lower branch of the 
levels pin to the value of El for high values of magnetic 
field B. At this point it would be useful to inspect the 
behavior of the two branches, E ni and E ri2 , Eq. (J7J for 
small and large values of magnetic field. For small mag- 
netic fields one has the following asymptotic expressions 



E n ± = -\E Cn + E L ± \j{E Cn -E L ) 2 + ixC% 



(7) 

The notation Ec n pertains to the pure Landau levels 
Ecn = (n+ j)huj c + h 2 k1/{2m*). From now on, index 



E m =-(E L -D) + A l (ni + -\ Ikv,. 



+Bi («i + §) (f^c) 2 + 0({huj c f) , (10a) 



E n 2 = ^(E L +D) + A 2 (n 2 + -) /;.,-. 
+B 2 (n 2 + 1 J (^ c ) 2 + 0{{huj c f) , (10b) 
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where I? = \/E\ -r i^^jy 

_Bi ; 2 are given by the expressions 



4:xCm and the coefficients Ai 2 and 
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These expressions are essentially the same as those of 
Ref.[la]. The overlap between the localized impurity 
wavefunction and conduction band states is neglected 
which is justified for B < B max = 35T. 



The reduction of the bandgap is manifested in Eq. (|10al) , 
as the zeroth order term, and is equal to (El — D)/2. The 
numerical values of A\ t2 , B\ 2 for the concentration of 
nitrogen impurities x = 0.08% are A\ = 0.92, A 2 = 0.08 
and B\ = —0.55 = — B 2 . On the other hand, we will also 
need approximate expressions for large B, 
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It is obvious that the lower branch E n2 pins to the 
nitrogen induced level, El for large magnetic fields. Fur- 
thermore, there is a gap between the values for B — » 
and B — > 00 which is given by 
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This is a consequence of the fact that we have a two-level 
like problem whose energy separation is determined by 
the matrix element Cjv, and the concentration x. For 
instance, for x — 0.08%, its value is So = 23meV. 

Raman measurements 14] indicate that the localization 
radius lies in the range cihw = 1-25 — 1.7nm (ajjw = 
a\/ln2), so that the parameter a is of order 2nm. For 
the present model to be valid, one should impose the 
condition that the cyclotron orbit is at least three times 
larger than the spatial extent of the impurity wavefunc- 
tion \i(jl\ 2 - This ensures that the magnetic field does not 
distort the impurity wavefunction, and it yields the upper 
bound B max = 35T for the magnetic field. The coeffi- 
cients a and /3 are found from their ratio (determined by 
the corresponding eigenvalue equation) and the condition 
that the wavefunction be normalized (the phase factors 
of a and /3 are taken to be zero) : 
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which results into 



III. ELECTRON-PHONON INTERACTION IN 

GAAS:N 

In order to obtain more precise values for the energy 
levels, one must also take into account the interaction 
of electrons with the thermal vibrations of the crystal. 
In a polar semiconductor like GaAs, electrons interact 
with longitudinal optical (LO) phonons more strongly 
than with other types of phonons. They may be assumed 
dispersionless having energy Elo — tv-^LO — 36meV in 
GaAs [lf|. Bearing in mind that coupling with the LO 
phonons is weak in common semiconductors, Frohlich[17J| 
proposed the following form of the Hamiltonian 



H = H + hujLO ■ J2 5 q 5 i + H < 



e—ph 



(17) 



where the expression for the interaction part with the LO 
phonons reads 



H e - P h = j2 v « eUir ( b * +b -J> 
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The strength of the electron-phonon interaction de- 
pends on the dimensionless coupling constant a which 
for GaAs has a small value of[16j 0.068. The interaction 
becomes important when the cyclotron energy Huj c ap- 
proaches the energy of the longitudinal optical phonon 
fvjjLO- This corresponds to the situation when, for in- 
stance, the energy of the unperturbed ground state with 



one re al p honon \n = 0, l p h), E { 



(0) 



hhJLO + 1/27jw c 



crosses [16J the first excited Landau state with no phonons 
\n = 1, Oph), E{ = (3/2)hu c . The electron-phonon 
interaction removes the degeneracy at uj c = lolo and 
makes the two level anticross. The renormalized values 
E , and E\ of the two levels can be found using perturba- 
tive methods. However, the usual Rayleigh-Schrodinger 
perturbation theory (RSPT) does not give precise values 
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FIG. 2: The first transition energy E\ — Eq in GaAs:N as 
a function of the cyclotron energy hw c , in units of the LO 
phonon energy ftwLO- It consists of two branches, one below 
the phonon continuum and another above it. For comparison 
the values of E\ — Eq for pure GaAs are shown by the dashed 
curves. 



The anticrossing behavior of the levels Eq and E\ is ob- 
vious. The lower curves, which correspond to the case 
below the phonon continuum pin to huiLO for large mag- 
netic fields. The difference between GaAs:N and GaAs 
is the largest when the cyclotron energy is comparable 
to the LO phonon energy, Hujlo- This difference should 
increase as the cyclotron orbit becomes comparable to 
the localization radius l c ~ o, but eventually both values 
should pin to hijjj^o for very large B. Another impor- 
tant difference between GaAs:N and pure GaAs is that 
the values of E\ — Eq above the phonon continuum in 
former (doped) case do not tend to hu c , but to a some- 
what smaller slope around A\ (see Eq. {TOa}). This is 
the consequence of the nonparabolicity of the conduc- 
tion band of GaAs:N. The minimum difference between 
the two levels is AE\§ = 5meV and it occurs at magnetic 
field B = 23T just above the longitudinal phonon energy 

frhJ LO - 



IV. CYCLOTRON RESONANCE ABSORPTION 



of the excited states for large values of the magnetic field. 
To overcome this problem, an improved Wigner-Brillouin 
perturbation theory[H [jjl [ig[ (IWBPT) is usually em- 
ployed to determine the pinning values of the renormal- 
ized levels. The energy correction for the n th state within 
this method is given by the expression [16j 
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and the matrix element |M nm | 2 has the following form 
for the first two energy levels 
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where s = (q±l c ) 2 /2. According to this perturbation 
scheme, the denominator D nm depends on the energy 
correction itself |18| . AE n , 



Cyclotron resonance measurements is a standard 
technique [20( to measure effective masses in bulk semi- 
conductors. In pure semiconductors, when scattering on 
defects and impurities can be neglected, and for the case 
of parabolic bands, absorption should ideally consist of 
a single sharp peak located at the cyclotron energy fiw c . 
However, in case of GaAs:N, it will be shown that the 
absorption linewidths are naturally broadened due to the 
change in the conduction band structure caused by the 
nitrogen impurities. Absorption is a measurable quan- 
tity which is determined by the oscillator strength that 
is defined byHMl 
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Si 



m*Efi 
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(23) 



In this work we will consider only transitions between ad- 
jacent levels, ni — 1 ~> ni, that are normally only possible 
transitions in cyclotron measurements. In the calculation 
of the matrix elements d ni>ni -i = {n\ — l|9 y |ni), only 
one term survives [15j when the total wavefunction, Eq. 
© is inserted in Eq. ([23]), and 
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(24) 



E m (-q z )-hw L o + &E n -AE , (22) 

so that the solution must be sought self-consistently. 
The correction to the ground state energy AEq is equal 
to its counterpart within Rayleigh-Schrodinger pertur- 
bation theory, i.e. AE = AE^ S . Fig. [5] shows the 
values of the first transition energy E\ — Eq as a func- 
tion of cyclotron energy hio c (i.e. magnetic field), by 
the solid curves. For comparison, the values of the same 
quantity are given for pure GaAs by the dashed curve. 



while for the states higher than the impurity state, 
El, the expression is the same with ni replaced by n<i. 
The coefficients a ni and a ni _i are given by the expres- 
sion Eq. (|15[) . In the next two figures, Fig. [3] and Fig. 0] 
we present the oscillator strengths for the first five tran- 
sitions between adjacent levels n — 1 — > n, within the 
lower and upper subbands, as a function of magnetic 
field B. The values of the oscillator strengths for the 
lower subbands decrease with increasing magnetic field B 
(Fig. [3]), due to the presence of coefficients a ni and a ni _i 
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FIG. 3: The oscillator strengths for the first 5 transitions 
between adjacent levels m — 1 — > m, below the impurity, 
localized level El as a function of magnetic field B. 



FIG. 5: Absorption coefficient a a t s in arbitrary units as a 
function of the incident energy E — tkj at room temperature 
T = 300K for two values of magnetic fields: B — 10T (solid 
curve), and B — 20T (dashed curve). 
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FIG. 4: The oscillator strengths for the first 5 transitions 
between adjacent levels n-2 — 1 — > ri2, above the impurity, 
localized level El as a function of magnetic field B. 



in Eq. (|2"4"|) . This can be explained by the fact that con- 
duction like states E ni acquire somewhat of a localized 
nature [15j as they approach the pinning value of E L . On 
the other hand, the strengths for upper subbands start 
from very small values indicating their highly localized 
nature at small magnetic fields. At the end, it would 
be useful to calculate the absorption coefficient a a ba as a 
function of the energy of the incident light, fku. Within 
the dipole approximation, the general formula reads [24] 



"..''■(/'-•I = — ^£ WM(/(iSi)-/(25/)) (25) 



n r 2m* 



i,f 



where n r is the refractive index, pj is the joint density of 
states, and haj is the energy of the incident light. The last 
factor in Eq. (|25l) is the difference in the Fermi-Dirac dis- 
tribution of the initial and the final state (Ei, Ef). Note 



that in Eq. (|2"5)l one should use the modified (joint) den- 
sity of states, pj, to take proper account of the presence 
of N impurities. The sum in the same equation will in 
practice terminate due to the finite Fermi level and tem- 
perature, while in the case of interest / = i + 1 due to 
selection rules. In the next figure, Fig. [5] the theoreti- 
cal estimate for the absorption coefficient for intraband 
transitions within the lower subband are shown as a func- 
tion of the incident energy E = fko at room tempera- 
ture T — 300K, for two different values of magnetic field: 
B = 10T (solid curve) and B = 20T (dashed curve). In 
both cases the electron concentration was kept constant 
at n c = 2 x 10 17 cm -3 , so that the position of the Fermi 
level changes with magnetic field. To achieve better clar- 
ity, a break point is introduced on the cc-axis between 16 
and 22meV. The broadening and asymmetry of the ab- 
sorption lines is a consequence of the band nonparabol- 
icity of E ni (k z ) introduced by the nitrogen impurities. 
For the same reason, the local maxima are not located 
at huj c but do depend on the quantum number m, since 
the equation E ni (k z ) — E ni -i(k z ) — hu> is not trivial. 



V. SUMMARY AND CONCLUSIONS 

In this work we presented a model to determine the en- 
ergy levels in dilute nitride GaAs:N under applied mag- 
netic field, taking into account electron-phonon interac- 
tion. The model is based on the Band Anticrossing Model 
for GaAs:N but here modified to include the interaction 
with the longitudinal optical phonons treated within sec- 
ond order perturbation theory. It is assumed that the 
impurity wavefunction has a Gaussian shape, spanning 
just a few lattice constants. The influence of nitrogen 
impurities are characterized by three parameters: local- 
ized level E L , the matrix element (7/v and concentration 
of the impurities x. The polaron correction to the energy 



levels results in a smaller transition energy than from 
pure GaAs. This difference is the largest around the LO 
phonon energy fuoLOi but eventually becomes zero for 
large magnetic fields. Furthermore it is shown that the 
oscillator strength for lower subbands (with respect to 
El) decrease with increasing of magnetic field as they 
approach the pinning value of El- The situation for the 
upper branch is the opposite, since they acquire some- 
what of an extended like nature. This should be revealed 
in the absorption measurement, in the sense that for only 
higher magnetic fields the intraband transition within the 
upper branch may contribute significantly to the absorp- 
tion spectrum. At the end, the theoretical estimates of 
the absorption coefficient are done for the lower subband 



at room temperature T — 300K, for two different values 
of magnetic field B. The absorption lines have asymmet- 
ric shape due to the non-parabolicity of the conduction 
band caused by the short range impurity potential. 
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